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Abstract 

Representations as determinants of M x M dimensional matrices are obtained for equal- 
time temperature correlators of the anisotropic Heisenberg XY chain. These representations 
are simple deformations of the answers for the isotropic XX0 chain. In the thermodynamic 
limit, the correlators are expressed in terms of the Fredholm determinants of linear integral 
operators. 

Introduction 

The description of the correlation functions for the models solved by means of Bethe ansatz 
is based on the representation for the correlators as the Fredholm determinants of linear 
integral operators. Such representations were obtained for the first time in EJ, for the sim- 
plest two point equal-time correlators for the one-dimensional impenetrable bosons model. 
Later they were generalized for the case of time-dependent correlators for the models being 
the free-fermion point of models solved by means of Bethe Ansatz (impenetrable bosons 
01 and the isotropic XX0 Heisenberg chain flEJ) and also for the case of finite interaction 
(|ij [?j for the one-dimensional Bose gas and ||for the XXZ chain). Such representations 
allow to write classical integrable equations for correlators which can be used, in particu- 
lar, to calculate the long time and large distance asymptotics for the correlation functions 

pjEj, Bg. 
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It was understood a time ago jH], |l^, Q that the language of the classical differential 
equations is natural for the description of the correlation functions of quantum integrable 
models. The recent progress in this direction is described in detail in the book for the 
example of impenetrable bosons based on the approach elaborated in [ p"o| , p"l| , |l2] , |l3} , |l4| |l5[ . 
The most important part of this approach is to consider the Fredholm determinant of the 
linear operator appearing in the representation of the correlators as a r-function for a new 
system of classical integrable equations (see also J9[). These linear operators should be of a 
very special form - so called "integrable" integral operators (see fj"i[ pof ). 

In this paper we obtain some determinant representations for equal-time temperature 
correlators for the anisotropic Heisenberg XY chain both for finite lattice and in thermo- 
dynamic limit. This model was introduced and studied for the first time in j^). Later it 
was investigated by many authors (see ^| and references therein). 

To compute the correlation functions we use a modification of the approach proposed in 
|25[ and based on using the integration over the Grassman variables and the corresponding 
coherent states. It is an essential point in our paper. We don't use, however, the functional 
integral as in . Using the coherent states simplifies the calculations leading automatically 
to the answers in the necessary form (see for example Appendix B where we reproduce the 
results [5) for the XXO chain using our method). 

The representations for the equal-time temperature correlators of the anisotropic XY 
chain as determinants of M x M matrices obtained in this paper being a direct generalization 
of the representations for the isotropic case [^, ^, ^] differ from them only by changing 
"Fermi" (or "Bose") weight in the kernels of integral operators for a weight depending on 
the anisotropy parameter. On the other hand these representations generalize the results 
for the anisotropic chain obtained in for zero temperature. The "integrability" of the 
integral operators appearing in the thermodynamic limit is evident. 

One should note that the isotropic XXO chain is the " free-fermion point" for the XXZ 
Heisenberg chain which is a model solved by means of the standard Bethe Ansatz. The 
anisotropic XY chain is the "free-fermion line" for the XYZ chain which is a model where 
the usual Bethe Ansatz doesn't work. Therefore it is very interesting from our point of 
view that the answers for the equal-time correlators have the form of a simple deformation 
of the answers for the isotropic case. We hope to give a corresponding description for the 
time-depending correlation functions in our next publication. 

Our paper is organized as follows. 

In the first section we describe the model and give the basic facts about the diagonaliza- 
tion of the Hamiltonian using the Bogolyubov transformation following the classical works 
I2H pjf . In section 2 we introduce the coherent states for isotropic XXO model and calculate 
the matrix elements of operators between these states. In section 3 we describe the coherent 
states for the anisotropic XY chain and consider their relations with the coherent states 
of the isotropic chain. In section 4 the simplest correlator is calculated. In section 5 we 
calculate the local spin correlator for the finite anisotropic lattice. In section 6 we obtain 
the results in the thermodynamic limit. In appendix A the basic facts about Grassmanian 
coherent states are given. In appendix B the derivation of the results for the isotropic case 
is presented. 



1 The XY Heisenberg chain 

The Hamiltonian of the XY spin chain describing the interaction between the nearest neigh- 
bors spins 1/2 placed in the sites of one-dimensional periodical lattice in a constant magnetic 
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field h is 

H = H + 7 Hi-ftS z , (1.1) 



where 

Ho 



1 

2 ^(°"mO-m + l +0'mO'm + l); (I- 2 ) 



2 

m— 1 

M 

2 



1 

Hi = -- ^(0"m0"m + l +CTmO"m + l); (1-3) 



and the third component of the total spin is 



M 



<■ ( L4 ) 



m— 1 

The total number of sites M is supposed to be even. Pauli matrices are defined as usual 

K, <&] = 2i6 mn e a ^al(a, /?, 7 = x, y, z); (1.5) 

zb ^- / x I ■ y \ 
a m = 2 1°"™ ± iCr mJ, 

with the periodical boundary conditions 

«T*r+i = «rf ■ (1-6) 

Due to the symmetries of the Hamiltonian the sign of the magnetic field (as the sign of 
the Hamiltonian) is not essential. We will assume that h > 0. 
The Jordan-Wigner transformation 

a m = exp{i7rQ(m- 1)}<t+; 

a+ = a m exp {i7rQ(m - 1)} , (1.7) 
introduces the canonical fermion fields a m , a+ on the lattice, 

[a+,a+] + = 0, [a m ,a+] + = S mn . (1.8) 
Operator Q(m) is the operator of number of particles on the first m sites of the lattice, 

m 
3 = 1 

where q m is the operator of number of particles in the site m: 

q m = a+a m = (J m a+ = i(l - o^). (1.10) 
The operator of the total number of particles, 

N = Q(M), (1.11) 

commutes with the operators Ho and S z but does not commute with the operator Hi and 
thus the total Hamiltonian H does not conserve the number of " a-fermions" . At the same 
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time the operator (— 1) N = exp{±i7rN}, anticommuting with the creation and annihilation 
operators 

[(-l) N ,a+}+ = [(-l) N ,a m } + = 0, (1.12) 
commutes with any bilinear in a m ,a+ operators, in particular, with the Hamiltonian 

[(-1) N ,H]=0. (1.13) 



Periodical boundary conditions (1.6) for the spins lead to the following conditions for 
the fermions: 

aM+i = (-l) N ai; ati+i = ai~(-l) N . (1.14) 



Introducing projectors P 



P ± = i(l±(-1) N ) 



(P ± ) 2 = P ± ; P + + P" = l; P + P~=P~P + = 0; (1.15) 
P ± a m = a m P*; [H,P ± ]=0, 
one can rewrite the Hamiltonian in the following form [^3| 

H = H + P + + H~P . (1.16) 

Both operators can be rewritten formally in the same form 

1 M 

H± = 2 X/ [( a ™ am+1 + a tn+i a m) + 7( a mai+i + o TO+ ia m )] + (1.17) 



2 

m— 1 



fl J CL^CLm 



m — 1 

T + 



2 

the only difference between H + and H~ being the boundary conditions: 
ajwr+i = — oi ; o^ +1 = — ap, for H + , 

OM+i = ai; a~l {+1 =a~l; for H~. (1-18) 

Hence the Fourier transformations to the momentum representation are different for 
these Hamiltonians. We denote the sets of permitted quasimomenta X + for H + and X~ 
forH": 

X ± = {p : exp{ipM} = T l, p6(-ir,7r]}, (1.19) 

or explicitly: 

X+ = { Pl =-«-^ + p, 1 = 1,2,...,*}, 

X- = {p, = -ir-^i, Z = 1,2,...,m}. (1.20) 

The corresponding formulae for the Fourier transformation can be written in the following 
form 

exp{-i7r/4} r-s 
a m = = > a p exp{i(m - l)p\ , 



P ex± 



+ exp{j7r/4} 

dm — 



M 

p6X± 



o+ exp{-i(m- l)p}, (1.21) 
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(the summation is taken over p £ X + for H + and over p £ X for H ) and 



exp{i7r/4} 



aim exp {-• i(m - l)p} . 



+ exp{— i7r/4} 



m— 1 

M 



M 



^ a+ exp{i(m - l)p} . 



The Hamiltonians H in the momenta representation are written as 



pex± 



£(p)a,pa p + ^^-(cLpat p + a-pdp) 



Mh 
2 ' 



(1.22) 



(1.23) 



where 

e{p) = h — cosp; r(p) = 7sinp. (1-24) 

Diagonalization of these Hamiltonians can be done using the Bogolyubov transformation 
(different for H + and H - ) leading to new canonical fermion operators A v and Ap 



A p = a(p)a p - P(p)a± p ; 
A+ = a(p)a+ + P(p)a-p, 



where 



a(p) = cos • 



Hp). 



Pip) = i 



and the angle 9(p) is defined by relations: 



cos%) = -g^; sin%) 



E(p) ■ 



E{p) = Ve 2 (p) + r 2 (p)>0, (p^0,7r). 



(1.25) 
(1.26) 

(1.27) 
(1.28) 



The momenta p = 0, 7r (appearing only in H ) should be treated separately. Following |2 
we put 

A = a - Aq = clq, 

An = ci-Tr ; yl^T = a^Jr, 



and 



£(0) = h - 1 = e(0) (E(Q) < for ft< 1) 
-B(tt) =h + l = e(vr) > 0. 

The Hamiltonians can be diagonalized as follows 

H± = E(p)A+A p + E±, 

pGX± 



(1.29) 
(1.30) 



where the "vacuum energy" is 



■I E 



(1.31) 



(1.32) 



p6X± 



(to calculate E one should take into account the definition (1.30)). One should note that 
for h < 1 the value Eq is not the ground state energy Eg , in this case Eg = Eq + e(0). 
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2 Coherent states for the XXO chain 



In this section we give some formulae for matrix elements of operators between the coherent 
states of the isotropic XXO chain which are necessary for the following calculations. The 
corresponding Hamiltonians H^xq (see ( |l.23| )) 

H xxo = X! e (P) ffl p a f ~ ~2~ ( 2 ' 1 ) 

p6X± 

will be denoted simply in thi s sect ion for the simplification. They are diagonal already 



in terms of the operators a p , ( 1.21 ). One should note, however, that the Fock vacuum 
(which is the same for H^), 

o m |0)=0, <0|a+=0 (m = 1,2,..., M), (2.2) 

a p |0)=0, <0|a+=0, (p € X ± ), (0|0) = 1, 

is the ground state only for h > 1. 

One introduces the coherent states (see Appendix A) different for the Hamiltonians H + 
and H~ 

\<f>,±) =expi \ l°>> ( 2 - 3 ) 



\±| = (0|expi rt a « \ ■ ( 2 - 4 ) 



q£X± 



The parameters <f> q ,<f)q (Grassman algebra elements) anticommute with other parameters 
and with all the operators a,a + . The main properties of the coherent states ( |2.3| ), are 
described in Appendix A. They are eigenstates for the operators a p and aj: 

a p \<t>,±) = cf> p \<j>,±) (p € X ± ), (2.5) 

(0*,±|a+ = 0;<0*,±| (PGX±). (2.6) 



The scalar product of the coherent states of one type is given by the usual formulae (A.4) 



,+\4>,+) =expi \ > ( 2 - 7 ) 



-) = expi <t>* q <t> q \ . (2.8) 



When it ca nnot ca use misunderstandings the sums on the right hand sides of the equations 
@, Q, (0), Q will be denoted a+0 = £ <4<t><i\ <f> + 4> = £ 4>p<t>v etc. 
The scalar products of the coherent states of different type are given by 

((/>*, +\ip, -) =exp < ^ ct)* p L pq ip q \ =exp{4>*Lip}, (2.9) 
\ P ex+, q £X- I 
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{*!>*, -It, +) =exp < ^ i]pL pq <t> q f=ex.p{il)*L<l>}, 

where the matrix element of the M x M matrix L are 

2 1 i t a - ii 

cot - — - - i ; 



M 1 - exp{-i(p - q)} M 



(2.10) 



It becomes evident if one takes into account that the Fock vacuum is the same for both 
types of state s and rewri tes t he scalar products in the "coordinate representation" using 
the formulae (1.21) and (1.22). For example, 



. -) = (0| exp I ^2 4>* m am > exp < ^ a+Vm > |0) 



where 



exp 



>cp{i7r/4} 



M 

PGX + 

exp{— in/ 4} 



^ </)p exp {~i(m - l)p} , 



M 



^ ip q exp {i(m — l)q} . 



qex- 



Consider now the matrix elements of the operator exp{aQ(m)} where (see (F9)) Q(m) 
is the number of particles operator at the first m sites of the lattice: 



(2.11) 



(2.12) 



Q = ^2 a t ai — a i<3pi,P2( m ) a P2 = a + Q(m)a, 



(2.13) 



(the quasimomenta pi and P2 here correspond to the states of the same type: pi,P2 £ X + 
or pi,p2 £ X~). Using the properties of the matrix Q(m) (evident in the coordinate 
representation) , 



Q (m) = Q{m); exp{aQ(m)} = I + (e a — l)Q(m), 
one has for the matrix elements between two states of the same type 

<0*,±jexp{aQ(m)}|^,±> = exp {<£*[/+ (e Q - l)Q(m)]0} 



(2.14) 



<?i*[5pg + (e Q - l)Q M (m)]</> 9 > , 



the matrix elements of the M x M matrix Q(m) being given as 

Qpq{m) = exp{-i(m - l)p/2}Qf q ' (m) exp{i(m - l)q/2}, 



1 sm 



m(p-q) 



}pq (m) , , 



M 



(2.15) 

(2.16) 
(2.17) 
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(for the diagonal matrix elements one should use the l'Hopital rule, Q pp (m) — Q P °J (m 



m/M). 



One has for the states of different types, analogously to (2.9) 



, ±| exp{aQ(m)}j0, =f) = exp {0* [L + (e a - l)Q(m)]0} = 
= exp< ^ 4>l[Lpq + {e a - l)Q Pq {m)]4> q 

I peX±,q£XT 



(2.18) 



where the matrix elements L pq and Q pq (m) are given by the formulae (2.10), (2.16) (one 
should note however that now if p £ X + then q £ X" and vice versa). In particular, one 
has for the operator exp{i7rQ(m)} entering the Jordan- Wigner transformation 

((?y,±|exp{i7rQ(m)}j0,=F} = exp {>*L(m)0} = 

= exp< ^ <t>* P L pq (m)<t) q \ , (2.19) 

where 

L pq (m) — exp{— imp}L pq exp{img}. (2.20) 
Turn now to the matrix elements ("form factors") of the local spins. Since 

P ± Om = °mP T (a = x,y, or a = ±), (2.21) 

we need only the matrix elements of th e loc al ope rator s <r^ between the states of different 
type. A direct calculation using ([Tt]), ( |l.22j ) and (2.18) gives : 



(</>*, ±|cr+ \tp, =f) = V'm(T) exp {<j)*L(m - 



= ^m(=F) exp < ^2 </> P L pq (m - l)ip q 

,±\<T^\i>,T) = cj>* m {±)exp{cj>*L{m - l)tp} . 
where we use natural notations 

exp{— in/4} 



<M±) 



?/> q exp{i(m- l)q}, 



exp{i7r/4} 



^ 0pexp{-i(m- l)p} . 



pGX± 



(2.22) 



(2.23) 



One should note that the matrices L (|2. lOh , L(m) (2.20) and Q(m) (2.16) are related by 



the following simple formula which can be proved by direct calculation: 

2J L piq (m)L qP2 = S P1P2 - 2Q P1P2 (m), (2.24) 



(pi,p 2 €X + , q€X or (pi,p 2 € X , q £ X ). 

Using these formulae one can reproduce the time- dep endent temperature correlation 
functions for the isotropic XXO chain obtained in [^, ^, |2(| . The corresponding calculation 
is given it in Appendix B. In the next sections we consider the anisotropic XY chain. 
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3 Coherent states for the anisotropic XY chain 



In this section the formulae for the matrix elements of operators between the coherent 
states of the anisotropic XY chain are given. We consider now two sets of c anoni cal fermion 
operators, a p ,a p and A Py A p , related by the Bogolyubov transformation (1.25). For each 
of these sets we introduce the coherent states ("old" and "new" ones) and calculate the 
scalar products of old and new states. 

Consider first the relations between the two vacuum states |0) and |0)) for the old and 
new sets 

a p |0>=0, <a+ = 0, Vp; (0|0) = 1, (3.1) 
A p \0))=0, ((A+ = 0, «0|0» = 1. (3.2) 
These states are related as follows: 



|o» = jv~ 1/a n + |o); «o| = N- 1/2 {o\n, 

where the operators Sl + and f2 are 

n+=exp{§y>(p)a+at p j; 
Q, = exp < i T~(p)a- P t 



and 



r (p) = tan 



Hp) 



(3.3) 



(3.4) 



(3.5) 



The normalization coefficient N can be represented as a determinant of a diagonal M x M 
matrix, 

N = (0|f2 + £7|0) = det(J + T), (3.6) 
where I is the identity matrix and 



T = diag(iT(p)). 



(3.7) 



Properties (3.2) for the states (3.3) can be easily checked using the commutation rela- 
tions 

[a p , fi + ] = T(p)£l + aL p and [fl, at] = r(p)a_ p fi. 

To calculate the normalization coefficient one makes use of the equation (A.6) and repre- 
senting 



N 



J d^e exp {-n + \t(p)(z;£ p + e-p^p)} 



Let us make the change of variables 

1 



(3.8) 



The Jacobian of this transformation is equal to 1, d£d£* = dcuduj* . One can easily check 
that 

e*£ = w"w; \^r( P ){C P C-p + xi. P i P ) = -Lo*Tu 1 (3.9) 
p 

where the matrix T is defined in ([T?]) . Hence, N = J dujdu* exp{-L>* (I +T)w} = det(/+T), 



which proves the relation (3j 
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Introduce now the coherent states for the new set of operators (compare with (A..2)) 



|X»=expj^A+X p j|0»; 

^ X;AP } 



«X*| = «0|exp J^X;A p y (3.10) 
The Grassman algebra elements X p , X* have the same properties as the old parameters 



CpiCp- The formulae ( A.2)-(A.10) are valid, of course, also for the new states. The direct 
calculation gives the following representations for the scalar products of old and new states 

<£* \X)) = N~ 1/2 exp J i r (p)(CC P - X-pXp) + a ''(P)CX P \ , 
I p p ) 

{{X*\0=N- 1/2 e X pl i^r(p)(e- P ep-^-p) + ^^ 1 (p)X;e P l. (3.11) 
I p p ) 

Considering the XY model it is necessary to introduce different sets of operators a p ,a p 
and A p , A^ corresponding to the sets of momenta X+ (for H + ) and X- (for H~), see 
(1.19), (1.20). Thus the new vacuums |fi}} will be also different for H + and H (unlike the 
old vacuum j0}). We will not usually mark in our notations this difference but one should 
have it in mind. 



Using the equation (3.11), one can calculate the matrix elements of the operators 
exp{ — f3H ± } diagonal in the new representations (different for H + and H~!) between 
old coherent states (also different for H + and H~); 

(f |e _i3H± |£) = Af" 1 e- /3J5 o ± det(7- JpT)e^ D ^, (3.12) 
where the variables ui p , uj* (p G X + for H + and p G X~ for H~) are defined by the equation 



(3.£), and diagonal M x M matrices Jp and D are 

= diag0>(p)); jp(p) = exp{-/3£(p)}, (3.13) 

D = diag(d(p)); d(p)= MM^pfl (3.14) 



(the matrix T is defined in (3.7)) 



To prove it we use the completeness of the new states 

(r|e _/3M± |0 = / dXdX*dYdY*x 



x(C\X))((Y*\e-^ \Y)){{X*\0exp{-Y*X - X*Y}, 
and note that the operators are diagonal in the corresponding new representations 



(1.31) 




({Y"\e- m \Y)) = e-^o exp { > e~^Y;Y p ) . (3.15) 
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After the integration over Y, Y* one gets 

«*|b- /,h± IO= / dXdx*{C\x)){{x*\i)x 



-PE± ) V"^ -BE(p)X*X„ 

x e exp < y e w » p 



The integral over X, X* can be ca lculated by means of the change of variables as in (3.8) 



Finally one gets the formula (3.12) 



In the following sections the results obtained here will be used to calculate the equal-time 
correlators. 



4 The simplest correlator for the XY chain 



In this section the partition function Z and the generating functional of the third compo- 
nents of local spins are calculated for the anisotropic chain. We begin by consi derin g the 
partition function. The initial representation is the same as in th isotropic case (B.l) 



Z==:(Z+ 



where 



Z% = Trexpj-Zffl*}; 
Z% = Trtexpl-^H+H-l) 



One gets the following re prese ntations for the contributions (see for example | 
matrix Jp is defined in (3.13)): 



Zf=e 



-^odet(I + J fj )= Y[ ( 2 cosh 



pGX± 



-PB. 



det(J - Jp) 



Yl I 2sin]j 



(4.1) 

(4.2) 
the M x M 

(4.3) 

(4.4) 



pex± 

To obtain, e.g., the fermionic contributions one should use the representation (A. 5) for the 
trace of operators: 



dYdY*((Y*\e 



-/3H a 



|y))exp{F*F}, 



and the representation ( p,15[ ) for the matrix element involved; after that one should calculate 
the Gaussian integral leading to the equality (4.3). The "bosonic" contrubutions can be 
calculated analogously. 

It is worth mentioning that it is sometimes convenient to represent the answer differently 
using the old coherent states, for example, 



z| = J d^r<rie _/3H± !e)ex P {rc}- 



By means of equation (3.12) for the corresponding matrix element one gets 
Zp = Ar _1 e" /3B ? det(7 - J p T)det(I + D), 



(4.5) 
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Z B = N~ 1 e~ l3E o det (I - JpT)det(I - D). 
Turn now to the simplest equal-time temperature correlator 

Z 



G{M) = ^Tr(e aQ<m) e-' 3H ) . 



(4.6) 



(4.7) 



As for the isotropic chain ((B^ 1,(1^9)) it can be represented as a sum of four contributions 

1 



G(m) 



2Z 



{Z F G p + Z F G f + Z B G B — Z B G £ 



where 



Z%G% =Trf e aQ(m) e - dH±N ) , 



Z±G± = Tr ( e 



One has for the contributions G^ B the following representations as determinants of 
M x M matrices 

G% = det (l+(e a - l)Q {0) (m)n F ) , (4.10) 

G% = det (/ - {e a - l)Q (0) (m)Q B ) . (4.11) 

Here the matrix elements of Q^°\m) are given by ( |2,16|) with p,q £ X + for G J B and 
p, g G X - for Gp B . Diagonal M x M matrices Qf and fis are given by the following 
formulae (p e X ± for Gp B ) 



aQ(ra) -/3H 3 



(-1)' 



(4.8) 



(4.9) 



Qf = D(J + D)- 1 = diag(w F (p)); 



wf(p) 



1-e 



ie(p) 



tanh 



(3E(p) 



n B = D{I - D)- 1 = diag(w s (p)); 



(the matrix £> is defined in ( 3.14 )). 

Explain now how to calculate the contribution G F : 



Z+G F = / d£d£*dr)dr)*{r)*\e 



ckQ(w) 



^H+|M e ^*f-{*r) 



(4.12) 



(4.13) 



(4.14) 



The matrix elements on the right hand side are given by (2.15) and (3.12). We change the 
variables 77* — > — rj*; since M is even, the integration measure is invariant. After that we 
change the variables as in (3.8) 



ft> = 7j(»7p + Pp = 7f (»7p + 



(4.15) 



As a result of this change of variables the measure remains invariant and 

rf (J + (e Q - l)Q(m)) V = p*{I + (e a - l)Q(m))p. (4.16) 



12 



The integration over u),tj* of the factors depending on these variables in (4.14) gives 
detD exp{— p*D~ 1 p}. One can calculate the remaining Gauss integral on p, p* taking into 
account the representation (vLa) for the partition functions and equation (4.12): 



Z+G F = Z+det(I + (e a - l)Q(m)fi F ). 



(4.17) 



One should note that the matrix Q(m) (2.16) differs only by an evident similarity trans- 
formation with a diagonal matrix from the matrix Q^(m). Thus the representation (4.1C) 
for is proved. The derivation of the contribution G F is almost the same, one should 
only use the momenta from the set X~ . To calculate the bosonic contributions one should 
use the property of the operator (— 1) N , 



(-irio = i-o. 



(4.18) 



in the repr esent ation similar to (4.14) which change evidently the calculations leading to 
the result ( J4 . 1 lj ) . 

We should make a remark about the equations (4.10), (4.11). Since the sets X + and 
X~ are symmetric under the change of momenta pi — > —pi (there is an exception, it is the 
momenta and n from the set X~; but the result is valid also in this case) one can rewrite 
the answer as (the sign "+" on the right hand side corresponds to Gf and the sign "-" 
corresponds to Gb) 



G% B = det(J± {e a - l)Q(m)n F , B ) 
l)f^Q(m)fi^) = 

- i)n^|Q(m)n^) 



= det(/ ± (e 

= det(/ ± (e a - l)n^lQ(m)n^) = ( 4 - 19 ) 

= det(J ± (e Q - l)Q(m)n F , B ). 

The bar means here the complex conjugation and u>f,b(—p) = u>f,b{p) 
To conclude we discuss some limiting cases. 

In the zero temperature limit (/3 — > oo) the "odd" contributions to the correlator are 
cancelled and one gets the representation for G(m) 



G{m) = det(J + (e Q - 1)Q (0) (m)fi ) (T = 0), 



where 



Oo = diag ( ojo (p) ) , coo (p) 



(4.20) 



(4.21) 



coinciding with the result obtained in [^7| . On the other hand for the isotropic case (7 = 0) 
taking into account that the angle 6(p) in the Bogolyubov transformation is 



0(p) = -wsignp, jpj < k F ; 
0(p) = O, |p|>fc F (7 = 0), 



(4.22) 



one gets for the weights 



w f (P) = - 



ws(p) 



coth 



1 — tanh 



(5e{p) 



Pe(p) 



(7 = 0). 



(4.23) 



Here fcf = arccos h is the Fermi momentum and e(p) is the dispersion of the XX0 chain 
(see (1.24)). Thus one reproduces the answers for the isotropic case ^, pif . 
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5 Equal-time correlators of the local spins 

Here we consider the equal-time correlation functions of the local spin operators (/3 = 1/T): 



G (ab) (m) = (<Tm + iCrJ>T = -x Tr((J^ +1 crJe dH ), a, b = +, 



(5.1) 



These correlators on a finite lattice of length M can be represented (as in (B.8)) as a sum 
of four contributions 



G 



(ab) 



1 

2Z 



7+ r~i(ab),+ . y- n (ab),- , y+ s-, 
ZjpKjp -\- &p\jp + ^5^5 



(ab), 



y- n (ab),- 
Z/ B Lr B 



(5.2) 



where 



± n (ab),± 



Z^G 



F F 



± n (ab),± 



= Tr a m+1 a-!_e 



b -/3H 3 



z bG b 



Tr 



(a b -/SH* / -innA 

I ^m+iO"ie (-1) I , 



(5.3) 



(partition functions Zp B are given by (4.3), (4.4)). For the contributions we obtain the 



following representations as determinants ofMxM matrices: 

G ( p +h± = G ( + _) ' ± = ^-det(J + Utt F + aCn F ) 
da 



(m > 0), 



G ^ +> ' ± = ^ tr " F = ^H WF(p) (m = 0) 



(5.4) 



G ( + _),± = 1- ^trfi F (m = 0), 



r ,(++),± _ r 



F 

( — ).± 



det(J + TO F -iaSfip) 



(m > 0), 



(m > 0), 



r .(-+),± 



-^trfis; G^+ ) ' ± = l + ^tr« s (m = 0), 



G 



( — ),± 



9q 



det(7 -t/fiir -iaSflp) 



(m > 0). 



(5.5) 
(5.6) 
(5.7) 
(5.8) 
(5.9) 



The diagonal M x M weight matrices Q,f and SIb are defined in ( [4,12] ), (4.13). Matrix 
elements of M x M matrices U, C, S are given by the following formulae 



U P1P2 = - exp | - p 2 )| Q^° 1 ) ,p 2 (r 



(5.10) 



(the definition of the matrix Q ( (m) is in (2.17)) 



C, 



1 m . . 
' ! " 2 = M cos "2 ^ +P2 ^' 



(5.11) 



(5.12) 



It is necessary to emphasize that the momenta numerating the matrix elements pi,P2 £ X^ 
for the contributions Gpg + an 
is true for the weight matrices. 



for the contributions Gp lb g + and pi,P2 £ for the contributions G^f^' i the same thing 
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Calculating the functions G^f^'* is red uced to ca lcula ting Gaussian integrals in the 



Grassmanian variables. For example, using (A.5) and (A.6) we represent 

d^ded V d V *( V '\a a m+1 a b 1 \ V )(C\e- ,3Ii+ \Oe^^ e ^ 



ZpG F 



Consider the calculation of the matrix element 



(5.13) 



Fab = {V* Wm + lVllv) = 
°m + l 



Let us use the formulae ( 2.22 ) for the matrix elements of spin operators and make the 
change of variables 



<€X+), 



fj* = e 2 p v *; fj p = e 2 



(5.14) 



One gets the representa tion for the matrix element using the new variables (we omit tildes 
over the new variables (5.14)): 



d 



where 



lu = rj*PLP( + CPLPr) - CC, 

IX 1 

l -(.Ql+<12) 



L e " 1 f i (Pl+P2) 



/+_ = CR+C, (R 



+ )<11<]2 



= — e 2 
M 



/++ = ivR~(; {R-) Pq = j^e~ 
f..=iCR-v r ; (R-) qP = -^e^-^ 



l -(p-q) 



(the bar means the complex conjugation). The diagonal matrix P has the form 

P = diag (e~^ p ^ . 

Use now the following identities 



imq /r\ c 1 \ ivm 

s = (2d m ,o - l)e ^ ; 



(5.15) 

(5.16) 
(5.17) 

(5.18) 

(5.19) 

(5.20) 

(5.21) 



^e imq L q 



( P eX + , q eX-, m = 0, l,...,M- 1), 



(PLPR+PLP) P1P2 = -(R+) P1P2 (1 ~ 2S m ,o) 
(R-PLP) P1P2 — {R-) P1P2 
(PLPR-) P1P2 — —(1 — 25 m ,o)(-R-) plP2 



(5.22) 

(5.23) 
(5.24) 
(5.25) 
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I t sho uld b e emp hasized that the matrix elements of the matrices on the left hand side of 
( fT23| ) and ( pT^ ) are 

(R+) qiq2 (qi,2 SX~) and (R-) pq (p 6 X+, q £ X"), 

and on th e right hand side (-R+)pij>2j C^-)piP2 (Pi,2 £ X + ). As a result, the matrix elements 
Fab (5.15) are represented as follows 







ex.p{ri*(I + U)l} + ar]*R+ri} 



F+- = (5 m , + —) exp{rj" (I + U)r] + ar/* R+r/} 
oat 



F ++ = —exp{f]*(I + U)r] + a>r)*S-ri} 
da 



F — = — exp{ri* (I + U)ri ~ a-q* S-r/} 
oat 



(5.26) 



Here the matrix elements of the M x M mat rix U are defined as U P1P2 = [PLP 2 LP) PlP2 — 
S P1P2 and this representation leads to (5.10) if one takes into account (2.24); the matrices 
R+, R+ are given by (5.17), (3.18) and (5.12); the matrix R+ is defined by the relation 

(R+) P1P2 = (R+) P1P2 (1 - 2S m<0 ), (5.27) 



(5.28) 



Now we put the expressions ( |5.12|) into (5^) using also ( 3.12 ). Then it is not difficult to 
get for the contribution G]7 + '' + 

z£G { f +) ' + = N~ 1 e-' 3E odet(I - J T)-^- J d£d£*di)drfx 

x exp{-?7*(J + U + aR+)ri + w* Dlu - rj*£ - £*rj}, (5.29) 



(besides (3.8), we changed the variabl es as in (5.14) and also we changed rj* — > —if)- 
Final ly we change the variables as in (4. IE), (i),lf) — > (p,p*); (£,£*) — > (ui,u)*). Using 
(4.16) and equalities 

11* (I + U + aR+)rj = p*(I + U + otC)p - ~(p*S-p* - pS-p), 



we integrate over oj,u>* and then over p, p* in (5.29). As a result we get the following 
representation as a determinant of a 2M x 2M matrix: 

Z+G ( ~ +h+ = N- 1 e~ l3E odet(I - J T)detDx 



x^det 1 / 2 
da 



aS- 



-B 



V B 1 



tS- J 
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where we denoted B = I + D 1 + U + aC. Using the generalized Gauss algorithm it is not 
difficult to check that 



/ 



det 



aS- 



-B 



-aSL / 



det 2 B + 0(a 



(it is valid also for m=0). It leads to the representation for G^T + ' 1 ' + . Analogously one can 
calculate other contributions. 

One should note that the expressions (5.4)-(j5JJ) for equal-time temperature correlators 



of the anisotropic XY chain differ from the corresponding expressions for the isotropic XXO 
chain only by changing the weights. For example, the corresponding modification for the 
"fermionic" weightis as before 



tanh 



(3e(p) 



iO{p) 



tanh 



PE(p) \ 



(5.30) 



6 The thermodynamic limit 

Now we consider the correlators for the anisotropic XY chain in the thermodynamic limit 
(the length of the chain M — > oo, the magneti c fie ld h is fixed). In this limit the partition 
functions Zp B are divergent and using (4.3), ( |4.4| ) one can estimate 



Zb/Z f < C A ; C = tanh 



< 1, 



where -E max is the maximal value of the quasiparticle energy E (p) (1.28). S o only the 
"fermionic" contributions survive in the thermodynami c lim it in ( 1.8) and (5.2). The de- 
terminants of the M x M matrices in the expressions (4.10), (5.4) for these contributions 
become in the thermodynamic limit the Fredholm determinants of the corresponding in- 
tegral operators. It is explained explicitly, for example, in p]. So we get the following 



The correlator G(m) (4.7) is given by 



G(m) = lim ^Tr (e aCi{m) e~ l3li ) = det (/ + (e a - 1)V) 



(6.1) 



At the right hand side, there is the Fredholm determinant of a linear integral operator 
acting on functions f(p) on the interval — n < p < n 



(Vf)(p)= / V(p,q)f(q)dq, 



(6.2) 



j mea ns the identity operator on the interval and the kernel V(p,q) is (see (4.1C), (2.17) 



V(p,q) 



1 sin 



2tt sin£^ 



(6.3) 
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The weight uJp{q) is equal to 



M,) = i(l-e^tanh«) 



(the angle 9(p) and the quasiparticle energy E( p) a re denned in ( 1.27 ), ( L2£ )). 
In the thermodynamic limit the correlators (5.1) 



G (ai,) (m) = lim l T r(^ +l( r ie -' 3H ) 



(6.4) 



(6.5) 



can be also repr esen t ed a s Fredholm determinants of linear integral operators on the interval 
[— 7r,7r]. From (5.4)-(5.6) one gets for m > 0: 



G 



(-+), 



(m) = G (+_) (m) = ^det(J - W + aC) 



d 



i:)a 
d 



G (++) (m) = G ( ~- ) (m) 
The kernels of the operators W, C and S are 

W(p,q) : 



det(7 - W - iaS) 



(m > 0), 
(m > 0). 



| i(p-q) Sin 

-e 2 — 



"l(p-g) 



-wf(?), 



(6.6) 
(6.7) 

(6.8) 

(6.9) 
(6.10) 

These answers differ only by changing the weight ( [5.30D from the answers for the isotropic 
XXO chain. 

One should note that the formal answers for the time-dependent correlation functions 
were obtained recently in |p9]| . In our next paper we hope to give more clear answers for 
the time-dependent correlators. 

We are grateful to the Russian Foundation of Fundamental Research (grants 95-01-00476 
and 96-01-00807) and INTAS (grant INTAS-RFFR 95-414). One of us (V.S.K.) thanks also 
TOO "Cyclone" for the support. 



C(p, q) = — cos tu F (q) 

N 1 • m (P + l) I \ 

S(p,q) = ^: sm ^ -UF{q)- 



Appendix A 

Consider a set of canonical fermion operators 

a q ,dq ([dp, a q ]+ = [dp, a+]+ = 0, [a p , a+]+ = S p>q ). 

We denote by M the number of operators a (or a + ) in the set. It is supposed that the Fock 
vacuum exists and has the following properties 



a q |0)=0; <0|a+=0, Vg; <0|0> = 1. 
We introduce the coherent states [BOl 



(A.l) 
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If,) =e X p|^a+e 3 | |0); 
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<ri = (c 9 *i = <o|cxp|^qa 9 | 



(A.2) 



where the summation is taken over the whole set. The parameters £ g , £j (Grassman algebra 
elements) anticommute among themselves and with all the operators a q ,a q . One should 
emphasize that the star in £* means only that the corresponding parameter is connected 
with a bra- vector; we don't consider involutions on the Grassman algebra; parameters £ and 
are entirely independent. The coherent states (A.2) are the eigenstates for the creation 
and annihilation operators, respectively 

a P io = &>iO; (e\4 = (e\e P - (a.3) 

One can easily calculate the scalar product of two coherent states using the commutation 
relations between a, a + , £, £*: 

(no = ex P { v e*£ q \ = ex P {ra (a.4) 



The trace of an operator O can be represented as an integral over the anticommuting 
variables of the matrix elements of the operator between the coherent states (see [pOl): 



TrO = J d£dCexp{t£}(r|0|a (A.5) 
and the expansion for the identity operator is given by 

1 = (-1) M / rf^expi-reilOn, (A.6) 



(we supposed that the number of sites M is even and the coefficient (— 1) in such formulae 
is usually omitted). We use the following notation 

1 1 
If an operator L has the following form 

L = S dp L pq a q = a + La, 

P,Q 

where L is a M x M matrix (with matrix elements L pq ) which can be diagonalized by an 
unitary matrix U, 

L = U + DU; U + U = UU + =I; D = diag(D q ), 
then the matrix elements of the operator exp L are given by 



<0|exp{L}|£) = exp<^ £ Cp(exp{L}) P9 \ = exp{£* expL£}, 
I p, i ) 



(A.7) 



and 

Trexp{L} = Trexp{a + La} = det[I + exp{L}], (A.8) 
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(there is the trace of the operator on the left hand side of the last formula and the deter- 
minant of the M x M matrix on the right hand side) . 

The last equality follows from a well-known formula for the Gaussian integral over the 
anticommuting variables 

J d£d£* exp{CK£r]*£, + £*r?} = expj-r?* K'^detK. (A.9) 

We need also to use a formula valid for an antisymmetric matrix A, 

Jdtexp{±tAt + nt} =exp|ir;4- 1 r 7 } VdrtA, A = —A T . (A.IO) 

Appendix B 

Here we give the derivation of the formulae for the time-dependent correlators for the 
isotropic XXO chain with finite number of sites M (M is supposed to be even). We hope 
that readers will appreciate the simplicity of the derivation of the results using the integra- 
tion over the Grassmanian variables (compared with paper |B|, ^). Below H denotes the 
Hamiltonian of the XXO model We begin with the calculation of the partition function Z. 
Following [ p3| , p5[ we represent it in the form 

Z = Trexp{-/3H} = i(Z+ + Zp + Z+ - Z B ), (B.l) 

where (5 = 1/T is the inverse temperature and the contributions on the right hand side are 
defined by the formulae 

Z% = Trexp-j-^H*} = e mh/2 det(I + J B ); 

Zg = Tr[(-1) N cxp{-/3H ± }]e' 3M ' l/2 det(J - J ), (B.2) 
where Jp is a diagonal matrix 

J fJ =dm g (e-^ p) ), (B.3) 

(we emphasize that p 6 X + for Zp B and p 6 X~ for Zp B ). To obtain these representations 
one should use the following relation 

(-l) N |0,±} = |-0,±>, (B.4) 

expf-zm^l^i) = e mh ' 2 \J^,±) = 



3MH/2 

e exp 



4e-^4> P ) |0>, (B.5) 



and rewrite, for example, Z^ as a Gaussian integral 

z ± = /mi/2 f d(f)d(f) * expl^*!/ - J p ](f>}, (B.6) 



(see (JAJ), (|A.5D). 

Consider now the simplest equal-time temperature correlator, 



G(m) = ^Tr [e aQ(m) e~ ^ , (B.7) 
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which can be represented in the following form 



where 



G(m) — —(Z F G F + Z F G F + Z B G B — Z B G B ), 



Z±G±=Tr( e aQ(m) e- flH± ), 
(-D N ) 



Ztd = Tr ( e 



aQ(ra) e -(3H a 



The contributions can be represented as determinants ofMxM matrices: 

G F = dot (/ + (e Q - 1)Q (0) (m)e F ) , 

G% = det (J - (e a - l)Q (0) (m)e s ) 
where diagonal matrices of the Fermi and Bose weights 9p and Qb are 

1 



(B.i 



(B.9) 



(B.10) 



Of = J/3 [I + Jfj] 1 = diag 
6s = Jp[I - J^] -1 = diag 



Qf(p) = 
Qb{ P ) = 



1 _|_ e /3e(p) 
1 



Z&e{p) - 1 



(B.ll) 



and the matrix Q^°\m) with matrix elements Qp^(m)was defined in (2.17) (p, q G X + for 
Gj B and p, q € X~ for Gp B ). To obtain these representations one should rewrite, for 
example, the contribution Z B G B as 



ZZGZ = e pMh/2 



>*V,-|e aQ(m) |-^->> 



use the expression (2.18) for the matrix element under the integral, calc ulate th e Ga ussian 
integral, perform the similarity transformation Q(m) — > Q^°'(m) (see (2.16), (2.17)) and 
extract the coefficient det(7 — Jp) from the determinant obtained. 

Now we turn to the time-dependent correlation function of the local spins; as usual 

Because of the translation invariance, the local spin correlators <Tm 2 (* 2 )> cr m 1 (ii) depend 
only on differences m = mi — mi, t = t% — t\: 



/^y(ab) i ,\ i itH a — itH b — /3H\ t . 

G (m, t) — — lr(e <r m+1 e a 1 e ), a, o = ±. 
Z 



(B.12) 



Due to the property (2.21) the correlator can be represented as a sum of four contributions: 
G {ah \m,t) = i [z+G^ b) - + + ZpG { ; b) <- + Z+G ( * b) < + - ^G^'") , (B.13) 

where 



± n (ab),± 



Z^G 



F^F 



± n (ab),± 



= Tr 



B^B 



= Tr 



( itH. ± a 
|6 Om + l 

( e cr„ 



itHT 6 -/3H d 



iHT b -/3H d 



(B.ll) 
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We begin with the correlator G' + \m,t). One can represent, e.g., the contribution 

Z+G { ~ +) = e f)Mh/2 ( dWdWilf'e*'*-*'** 



x (0* J(0, t), +k" +1 i J(0, — *)^r, -}<^*, J fl 0, +), 
where we introduced a diagonal M x M matrix J(m,t): 

J(m, t) = diag ( e -*p™+"(p )*) . (B.15) 



Now we use an expression (2.22) for the matrix elements and perform the integration over 
xp, ip* and then over <f>, <j>*: 



da 

= e /3M,l/2 -|^det [/ + L(m, t)L(Q)J p + aRJp] 
where we defined the matrix 

L(m, t) = J(0, t)L(m) J(0, -t) 

(in the matrix elements L Pt q(m,t), p G X + and g G -X" - ), an d also the matrix R of rank 1 
(all the columns of this matrix are equal): 

R - _L„-»Pi m +^(pi)t p v-+ 

"pi ,y 2 — TV^" ' ' ' 



After the similarity transformation with the diagonal matrix J(— — |) one gets 

(B.16) 



Z+G<T +) ' + = e^ 72 A det [/ + LJ^ + a £j 5 l 

0Q L J 



where matrix elements of matrices L and R are 

L P iP2 = e lmpi/2 - Ite(pi)/2 ^L pl<? (m,t)L 9P2 (0)e- lmp2/2+,te(p2)/2 , (B.17) 

^ _ _j_g-!'rip 1 /2+iic(pi)/2 e -imp 2 /2+ite(p2)/2 jg-j 

(here pi,p2 £ X + ,q £ X~). 

To represent these matrices in a more convenient way we introduce the following func- 
tions (we follow the paper j^], changing a bit some notations; recall that e(q) = h — cosq) 



/ , \ 1 X ^ imq4-it cos a /t> i r\\ 

i imq-\-it cos g 



9 2 



22 



d{m,t,p) = jjJ2 



,2 q-p 



(B.21) 



Using identities 



n piP2 = e a{m, t,pi)eb(m, t,p 2 ) ((a, b) = +, -). 

^-^W = M 2 ( q ex-,pex + y, 

' — ' sin -2^- 



e_ (m, t,p) = e - ! ^ p - l 3 » s " ; (B.22) 

e^m,t,p) = e~(m,t,p)e(m,t,p). (B.23) 

All these functions are of order O(l) as M — > oo. It is also convenient to introduce four 
one-dimensional projectors IT ++ , II" 1 , II _ + , II which are M x M matrices with matrix 
elements 

(B.24) 

(B.25) 
(B.26) 

(B.27) 

(B.28) 
(B.29) 

(B.30) 



cot(:r — u) cot(:r — v) = cot(tt — u)[cot(a; — u) — cot(x — v) 
one can represent the matrix L in the following form 



L = i+ i- [5 + J n+-- in - + ] 



M 

The diagonal and nondiagonal matrix elements of the matrix S are given by 

S pp = d(p)e- lmp - ltcosp 



S » 



The matrix R is a projector, 



e+(pi)e-(p 2 ) - e-(pi)e+(p 2 ) 



tan^^a 



M 



(For brevity, the arguments m and t are omitted, being the same for all the functions 
d,e,e±). 

Other contributions can be calculated analogously. We give only the results: 



G ( - +) ' ± = e lth -^det 



d_ 

da 



I + —{S + iU + - - iU- + + an~~)@ F 



(B.31) 



G (-+),± = e itu 

aa 



i + —(s + m+- - m-+ + an— )6s 



Though formally both contributions G 



(.-+), ± 



(as G F + ' ,± ) are written in the same form, 



really they are, of course, different. It is necessary to take into account that 
p,p 1 ,p 2 eX + , q g X~ for G<T+ ) ' + ; 

p,Pi,p 2 €X~, qeX + for G ( ~P'~; 

in all the formulae ( B.1S| )-(B.24), (B.28) and ( B,2S| ) defining the functions g,e,d,e± and 
the matrix elements of the matrices S and II ab . 

Analogously one can calculate also the correlator G^ K The corresponding contribu- 
tions are 



n (+-),± _ ith 



9(m,t) + - 



dot 



23 



- a{n ++ + igU + ~ - igW + + g 2 n~})0 F 



(B.32) 



G 



(+-),± 



det 



*{n ++ +ig!l + ~ -igW 



h + 3 2 n— })e B 



(one should take into account once again the difference between the contributions G F i ~ B * > ' + 
and G^pg^'~ explained above). 

One should underline that the results obtained here coincide with the representations 
obtained using another method (for which rather complicated calculations are required) 
j] . Precisely our e xpres si on ca n be obtained by choosing arbitrary constants C\ = 
-i in the formulae (A.l)-( A..10) of the Appendix A of the paper |2SJ]. 
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